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Abstract 

> 

We present a generalization of the resonant neutrino conversion in matter, including a 
random component in the matter density profile. The study is focused on the effect of such 
matter perturbations upon both large and small mixing angle MSW solutions to the solar 
neutrino problem. This is carried out both for the active-active v e — > z/ MjT as well as active- 
sterile v e —> v s conversion channels. We find that the small mixing MSW solution is much more 
stable (especially in 5m 2 ) than the large mixing solution. Future solar neutrino experiments, 

> 

X 
S3 



such as Borexino, could probe solar matter density noise at the few percent level. 
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1. The comparison among the present experimental results on the observation of the solar 
neutrinos strongly points to a deficit of neutrino flux (dubbed the Solar Neutrino Problem 
(SNP)). The most recent averaged data of the chlorine 0, gallium || |J and Kamiokande 
experiments are: 

R%* = (2.55 ± 0.25)SNU, R^l = (74 ± 8)SNU 1 , R^l = (0.44 ± 0.06)i?ff 5 (1) 

where R^a 95 * s the prediction according to the most recent Standard Solar Model (SSM)by 
Bahcall-Pinsonneault (BP95)§ . 

In particular the SNP is now understood as the strong deficit of the beryllium neutrinos 
J7J. On the other hand, the high energy boron neutrinos are moderately suppressed, while 
the low energy ones are almost undepleted. All this seems to imply that any astrophysical 
solution fails j7|, [| in reconciling the experimental data. 

From the particle physics point of view, however, the resonant neutrino conversion due 
to the neutrino interactions with constituents of the solar material (the Mikheyev-Smirnov- 
Wolfenstein (MSW) effect) Q offers the best explanation of the present experimental situation. 
This scenario provides an extremely good data fit in the small mixing region with 5m 2 ~ 



10" 5 eV 2 and sin 2 20 ~ KT 3 1(T 2 |I|P,||. Moreover, the study of the MSW effect has 



revealed its stability against possible changes of the SSM input parameters |TT| especially in 
the Sm 2 parameter. 

This talk deals with the stability of the MSW solution with respect to the possible presence 
of random perturbations in the solar matter density M . 



We remind that in Ref. |H| the effect of periodic matter density perturbations added to a 
mean matter density p upon resonant neutrino conversion was investigated. The major effects 
show up when the fixed frequency of the perturbation is close to the neutrino oscillation eigen- 
frequency, and for rather large amplitude values (~ 10 -r- 20%), giving rise to the parametric 
effects ||13|| . There are also a number of papers which address similar effects by different 



approaches [TJ, ITS] . 

Here we consider the effect of random matter density perturbations Sp(r), characterised 
by an arbitrary wave number k, 

Sp(r) = J dkSp(k) sin kr , (2) 

rather than a periodic or regular perturbation. The effect of solar density as well as solar 
magnetic field fluctuations upon neutrino spin-flavour conversions has also been considered in 



Ref. [15 1, using somewhat different methods. 



Moreover, as in Ref. ||15||, we assume that the perturbation 8p has Gaussian distribution. 
For small inhomogeneities, the spatial correlation function (£ 2 ) can be taken as 



(Sp(r 1 )Sp(r 2 )) = 2p 2 (e)L S(r 1 -r 2 ), (£ 2 > = M , ( 3 ) 

P 



whose correlation length L obeys the following relation: 



I free *C Lq ^ X m 



(4) 



where Zf ree = (crno)" 1 is the mean free path of the electrons in the solar medium^, and X m 
is the neutrino matter wave length. The lower bound, is dictated by the hydro-dynamical 
approximation used later, whereas the upper bound expresses the fact that the scale of fluctu- 
ations should be much smaller than X m (the characteristic neutrino propagation length scale), 
as indeed the Eq. (|3|) requires. For the sake of discussion, in the following we choose to adjust 
L as follows: 

io = 0.1 x A m . (5) 

The SSM in itself cannot account for the existence of density perturbations, since it is based 
on hydrostatic evolution equations. On the other hand, the present helioseismology observa- 



tions cannot exclude the existence of few percent level of matter density fluctuations [16], [17 . 
Therefore, in what follow we assume, on phenomenological grounds, such levels for £, up to 
8%. 

Before generalizing the MSW scenario, accounting for the presence in the interior of the 
sun of such matter density fluctuations, first we give a quick reminder to the main features of 
the MSW effect. 

2. The resonant conversion of neutrinos in a matter background is due to the coherent 
neutrino scattering off matter constituents |J. This determines an effective matter potential 
V for neutrinos. In the rest frame of the unpolarised matter, the potential is given, in the 
framework of the Standard Model, by 

V = (6) 

m p 

where Gf is the Fermi constant and Y is a number which depends on the neutrino type and on 
the chemical content of the medium. More precisely, Y = Y e — ^Y n for the v e state, Y = —\Y n 
for Up and v T and Y = for the sterile v s state, where Y e ^ n denotes the electron and neutron 
number per nucleon. Let us note the dependence of V on the matter density p for which one 
usually consider the smooth distribution, as given by the SSM || [L8|, [19| . 

Once assumed that there exists a non-vanishing mass difference 5m 2 between two different 
neutrino states and a non- vanishing neutrino mixing 9 in vacuum, the neutrinos z/ e 's, created in 
the inner region of the sun, where the p distribution is maximal, can be completely converted 
into u y (y = p, t or s), while travelling to the solar surface. 
This requires two conditions ||: 

1) - the resonance condition. Neutrinos of given energy E experience the resonance if the 
energy splitting in the vacuum 5m? cos 26 /2E is compensated by the effective matter potential 
difference AV ey = V e — V y . It is helpful to define the following dynamical factor A ey 

A ey (r) = ±lAV ey (r)- 5 -^ C os26] (7) 

which vanishes at the resonance, A ey = 0. This condition determines the value 

pres — ( m p cos 29/2 a/2Gf) (Y" e — Y y )5m 2 / E which, in turn, implies a resonance layer Ar. 

2 For Coulomb interactions, the cross-section a is determined by the classical radius of 
electron ro e = e 2 /m e c 2 ~2x 10~ 13 cm, resulting in Zf rco ~ 10 cm for a solar mean density 
n ~ 10 24 cm" 3 . 



2) - The adiabatic condition. At the resonance layer, the neutrino conversion v e — > v y 
is efficient if the propagation is adiabatic. This can be nicely expressed requiring that the 
neutrino wavelenght A m be smaller than Ar 0, 

* //> x 5m 2 sin 2 29R Q „ „ „ „ /„x 

a r = Ar/(A m ) res = — — — ^ > 1 , i? o ~O.LR , (8) 

47T£/ COS Z(7 

7T 

A m = ; = , Ar = 2p res tan26'|dp/dr|~ 1 . 

lA 2 ey + {5m 2 ) 2 sin 2 29/16E 2 



2. Now we re-formulate the evolution eqaution for the neutrino accounting for a fluctu- 
ation term Sp superimposed to the main profile p. The perturbation level £ = — induces a 
corresponding random component for the matter potential of the form AV ey ^. The evolution 
for the u e — v y system is governed by 

d_ (v e \ ( H e H ey \( v e \ 
% dt \v y ) \ H ey H y ) \v y ) ' W 

where the entries of the Hamiltonian matrix are given by 

~ 5tyi 2 
H e = 2[A ey (t) + A ey (t)\, H y = 0, = — sin 29, 

I Stti 2 ~ i 

A ey (t) = ~[AV ey (t) - — cos20], A ey (t) = ~AV ey {t)t (10) 

Here the matter potential for the active-active v e — > v^ T conversion reads 

AV eKT) (t) = ^^p(t)(l-Y n ) (11) 



m p 



or alternatively in case v e 



AV es (t) = ^^p(t)(l-lY n ) (12) 



(the neutral matter relation Y e = 1 — Y n has been used). 

The above system can be rewritten in terms of the following equations: 

P(t) = 2H ey I(t) 
R(t) = -H e (t)I(t) 

I(t) = H e (t)R(t) - H ey (2P(t) - 1) (13) 

where P = \v e \ 2 is the v e survival probability, R = Re(z/ y z/*) and I = lm.(v y v*) with the 
corresponding initial conditions P(t ) = 1, I(t ) = 0, R(t ) = 0. The Eqs. (|13|) have to be 
averaged (see |L| for more details) over the random density distribution, taking into account 
that for the random component we have: 



= 0, (A ey {t)A ey {t x )) = K8(t - 10, (14) 
1 

— L 

2 



1 



«(*) = (K(t)) L o = o^ y (t)(e)L . (15) 



The noise-averaged version of the system flUf ) reads as 



V(t) = 2H ey X{YA) 

K(t) = -2A ey (t)l(U) - en(U)n(U) 

±{t) = 2A ey (t)TZ(U) - e/c(U)J(U) - H M {eV{u) - oo). (16) 

where clearly (P(t)) = V(U), (R(t)) = TZ(U), (I(t)) = T(U). As expected the system of 
equations (|i~6"D explicitly exhibits the noise parameter k. 

It is now possible to envisage the main effects due to the presence of the random field 5p 
upon the MSW scenario. Now the " dynamics " is governed by one more quantity i.e. the 
noise parameter k, besides the factor A ey . Actually, the quantity k can be given the meaning 
of energy quantum associated with the matter density perturbation. However, let us note that 
the MSW resonance condition, i.e. A ey (t) = remains unchanged, due to the random nature 
of the matter perturbations. The comparison between the noise parameter k in Eq. (|i"5|) and 
A ey (t) shows that /t(t) < A ey (t), for £ < few %, except at the resonance region. As a result, 
the density perturbation can have its maximal effect just at the resonance. Furthermore, 
one can find the analogous of condition 2) (see Eq. (^) for the noise to give rise to sizeable 



effects. Since the noise term gives rise to a damping term in the system (|16|), it follows that 
the corresponding noise length scale 1/k be much smaller than the thickness of the resonance 
layer Ar. In other words, the following adiabaticity condition 

a r = Ar n res > 1, (17) 

is also necessary. There is a simple relation between the two adiabaticity parameters a r (cfr. 
(§)) and a r : 

& ^ ar ^kw- (18) 

For the range of parameters we are considering, £ ~ 10 -2 and tan 2 26 > 10~ 3 — 10~ 2 , and due 
to the r.h.s of @, there results a r < a r . This relation can be rewritten as K res < 5H res , where 
SH res is the level splitting between the energies of the neutrino mass eigenstates at resonance. 
This shows that the noise energy quantum is unable to "excite" the system, causing the level 



crossing (even at the resonance) [13]. In other words, it never violates the MSW adiabaticity 
condition. From Eq. ( |IBD it follows also that, in the adiabatic regime a r > 1, the smaller the 
mixing angle value the larger the effect of the noise. Finally, as already noted above, the MSW 
non-adiabaticity a r < 1 is always transmitted to a r < 1. As a result, under our assumptions 
the fluctuations are expected to be ineffective in the non-adiabatic MSW regime. 

3. All this preliminary discussion is illustrated in the Fig. 1. For definiteness we take 
BP95 SSM H as reference model. We plot V as a function of E/5m 2 for different values 
of the noise parameter £. For comparison, the standard MSW case £ = is also shown 
(lower solid curve). One can see that in both cases of small and large mixing (Fig. la and 
Fig. lb, respectively), the effect of the matter density noise is to raise the bottom of the pit 
(see dotted and dashed curves). For example, the enhancement of the survival probability 
can easily reach 20% for £ values as small as 4%. In other words, the noise weakens the 
MSW suppression in the adiabatic-resonant regime, whereas its effect is negligible in the non- 



adiabatic region, in complete agreement with the results of Ref . |L5f . The relative increase of 



the survival probability V is larger for the case of small mixing (Fig. la) as already guessed on 
the basis of Eq. (|l"8"l). We have also drawn pictorially (solid vertical line) the position, in the 
V profile, where 7 Be neutrinos fall in for the relevant 5m 2 ~ 10~ 5 eV 2 , to visualize that these 
intermediate energy neutrinos are the ones most likely to be affected by the matter noise. 

4. Let us analyse the possible impact of this scenario in the determination of solar neutrino 
parameters from the experimental data. For that we have performed the standard x 2 fit in 
the (sin 2 29, 5m 2 ) parameter space. The results of the fitting are shown in Fig. 2 where the 
90% confidence level (C.L.) areas are drawn for different values of £. Fig. 2a and Fig. 2b 
refer to the cases of v e — > v^ T and v e — > v s conversion, respectively. One can observe that 
the small-mixing region is almost stable, with a slight shift down of 5m 2 values and a slight 
shift of sin 2 29 towards larger values. The large mixing area is also pretty stable, exhibiting 
the tendency to shift to smaller 5m 2 and sin 2 29. The smaller 5m 2 values compensate for 
the weakening of the MSW suppression due to the presence of matter noise, so that a larger 
portion of the neutrino energy spectrum can be converted. The £ = 8% case, considered for 
the sake of demonstration, clearly shows that the small mixing region is much more stable 
than the large mixing one even for such large value of the noise. Moreover the strong selective 
7 Be neutrino suppression, which is the nice feature of the MSW effect, is somewhat degraded 
by the presence of matter noise. Consequently the longstanding conflict between chlorine and 
Kamiokande data is exacerbated and the data fit gets worse. Indeed, the presence of the 
matter density noise makes the data fit a little poorer: Xmin = 0.1 for £ = 0, it becomes 
Xmin = 0-8 f° r £ = 4% and even Xmin = 2 for £ =8% for the v e — > v^ T transition. 

The same tendency is met in the case of transition into a sterile state (Fig. 2b): Xmin = 1 
for £ = 0, it becomes Xmin = 3-6 for £ = 4% and Xmin = 9 for £ =8%. 

In conclusion we have shown that the MSW solution to the SNP exists for any realistic 
levels of matter density noise (£ < 4%). Moreover the MSW solution is essentially stable in 
mass (4 ■ 10" 6 eV 2 < 5m 2 < 10" 5 eV 2 at 90% CL), whereas the mixing appears more sensitive 
to the level of fluctuations. 



5. We can reverse our point of view, wondering whether the solar neutrino experiments 
can be a tool to get information on the the level of matter noise in the sun. In particular, 
the future Borexino experiment ]20|], aiming to detect the 7 Be neutrino flux could be sensitive 
to the presence of solar matter fluctuations. In the relevant MSW parameter region for the 
noiseless case, the Borexino signal cannot be definitely predicted (see Fig. 3a). Within the 
present allowed C.L. regions (dotted line) the expected rate, Z Be = R P Be d / Rse 95 (solid lines), 
is in the range 0.2 0.7. 

On the other hand, when the matter density noise is switched on, e.g. £ = 4% (see 
Fig. 3b), the minimal allowed value for Z^ e becomes higher, Z^ e >0.4. Hence, if the MSW 
mechanism is responsable for the solar neutrino deficit and Borexino experiment detects a 
low signal, say Ze e < 0.3 (with good accuracy) this will imply that a 4% level of matter 
fluctuations in the central region of the sun is unlikely . The same argument can be applied 
to v e — > v s resonant conversion, whenever future large detectors such as Super-Kamiokande 
and/or the Sudbury Neutrino Observatory (SNO) establish through, e.g. the measurement of 
the charged to neutral current ratio, that the deficit of solar neutrinos is due to this kind of 



transition. The expected signal in Borexino is very small ~ 0.02 for £ = (see Fig. 3c). 
On the other hand with £ = 4%, the minimum expected Borexino signal is 10 times higher 
than in the noiseless case, so that if Borexino detects a rate Z Be < 0.1 (see Fig. 3d) this would 
again exclude noise levels above 4%. 

Let us notice that Super-Kamiokande and SNO experiments, being sensitive only to the 
higher energy Boron neutrinos, probably do not offer similar possibility to probe such matter 
fluctuations in the sun. 

The previous discussion, which certainly deserves a more accurate analysis involving also 
the theoretical uncertainties in the 7 Be neutrino flux, shows the close link between neutrino 
physics and solar physics. 

It is a pleasure to thank N. Yahlali and H. Nunokawa for reading the manuscript. This 
work has been supported by the grant N. ERBCHB1 CT-941592 of the Human Capital and 
Mobility Program. 
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Fig. 1: The averaged solar neutrino survival probability P versus E/Sm 2 for small mixing an- 
gle, sin 2 29 = 0.01, (Fig. la) and for large mixing angle, sin 2 29 = 0.7, (Fig. lb). The different 
curves refer to different values of matter noise level £ as indicated. 
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Fig. 2: The 90% C.L. allowed regions for the v e — > z/ MjT (Fig. 2a) and for the z/ e — > z/ s (Fig. 2b) 
conversion. The different curves refer to different values of matter noise level £ as indicated. 
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Fig. 3: The iso Z Be = R P bI /-Rff 95 contours (figures at curve) in the v — e scattering Borexino 
detector (solid lines). The threshold energy for the recoil electron detection is 0.25 MeV. The 
90% C.L. regions (dotted line) and the corresponding best fit point are also drawn. Fig. 3a 
and Fig. 3b refer to the case of v e — > v^ T conversion and for £ = and £ = 4%, respectively. 
Fig. 3c and Fig. 3d refer to the case of v e — > v s conversion and for £ = and £ = 4%, 
respectively. 



